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We use a general circulation model (NASA Goddard Institute for Space Studies GCM 3) to investigate the 
regional climate response to removal of aerosols over the United States. We perform a pair of transient 
2010-2050 climate simulations following a scenario of increasing greenhouse gas concentrations, with 
and without aerosols over the United States and with present-day aerosols elsewhere. We find that 
removing U.S. aerosol significantly enhances the warming from greenhouse gases in a spatial pattern that 
strongly correlates with that of the aerosol. Warming is nearly negligible outside the United States, but 
annual mean surface temperatures increase by 0.4-0.6 K in the eastern United States. Temperatures 
during summer heat waves in the Northeast rise by as much as 1—2 K due to aerosol removal, driven in 
part by positive feedbacks involving soil moisture and low cloud cover. Reducing U.S. aerosol sources to 
achieve air quality objectives could thus have significant unintended regional warming consequences. 

© 2011 Elsevier Ltd. All rights reserved. 


1. Introduction 

Atmospheric aerosols affect climate directly, through scattering 
and absorption of sunlight, and indirectly, by providing cloud 
condensation nuclei. Increasing anthropogenic aerosols have 
exerted a significant cooling effect over the past century (Forster 
et al., 2007; Myhre, 2009; Murphy et al., 2009), mitigating in part 
the warming due to greenhouse gases (e.g., Hansen et al., 2007). 
U.S. emissions of aerosols and their precursors are presently 
declining sharply due to policies aimed at air quality and acid rain 
(Streets et al., 2004; Environmental Protection Agency, 2006). The 
regional climate effect could potentially be large since the present- 
day direct radiative forcing over the United States from anthropo- 
genic aerosols is about -2 W nrT 2 , similar in magnitude but 
opposite in sign to the globally averaged anthropogenic greenhouse 
forcing (Liao et al., 2004; Koch et al., 2007; Shindell et al., 2007; 
Forster et al., 2007). However, the connection from regional 
forcing to regional climate response is not well understood 
(National Research Council, 2005). We use here a general circula- 
tion model (GCM) in a 2010-2050 transient-climate simulation to 
investigate the climate response to suppression of U.S. aerosol 
sources. 
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Ground-based observations of clear-sky solar radiation over the 
United States reveal a surface dimming of -2 W nrr 2 per decade 
over 1960—1990 (Liepert and Tegen, 2002). This trend reversed to 
a brightening of +5 W m~ 2 per decade between 1995 and 2007 
(Long et al., 2009). U.S. trends in the sources of anthropogenic 
aerosols roughly follow these trends in insolation (Lefohn et al., 
1999; van Aardenne et al., 2001; Streets et al., 2006; Bond et al., 
2007). For example, the EDGAR inventory shows U.S. sulfur emis- 
sions increasing by 30% from 1950 to 1980 and then returning 
tol950s values over the next two decades (van Aardenne et al., 
2001; Olivier and Berdowski, 2001). Volcanic aerosol from the 
eruptions of El Chichon (1982) and Mount Pinatubo (1991) likely 
contributed to transient surface cooling but with little impact on 
the above long-term trends (Wild et al., 2005; Wild, 2009). 

Observations suggest that declining aerosols over Europe 
contributed to an observed continental surface warming of 1 K for 
1980-2006 (Ruckstuhl et al., 2008; Philipona et al., 2009). All-sky 
solar radiation at the surface increased by about +2 W nrT 2 
during this time period in Europe, concurrent with a 60% decline in 
regional aerosol optical depths (AOD) (Ruckstuhl et al., 2008). 
Philipona et al. (2009) found that the 1980-2006 shortwave 
forcing due to aerosol decline was 2-3 times larger over the 
continent than the concurrent longwave forcing due to rising 
anthropogenic greenhouse gases. Over India, the aerosol burden 
has increased dramatically in recent decades, yielding surface 
forcings as large as -15 to -25 W nrr 2 in some regions during the 
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dry season (Krishnan and Ramanathan, 2002; Ramana et al., 2004). 
Krishnan and Ramanathan (2002) provide observational evidence 
of a 0.5 K cooling over India since the 1970s due to increasing 
aerosols. Over the United States, surface temperatures rose -1 K 
from 1980 to 2006 after having remained flat from 1960 to 1980 
[http://www.ncdc.noaa.gov], though no study to date has attemp- 
ted to relate this rate of warming to U.S. aerosol forcings. 

In a closed system, regional radiative forcing from anthropogenic 
aerosols would cause a collocated climate response over the source 
continent. Horizontal transport of heat dilutes the response (National 
Research Council, 2005). Several GCM studies have reported that the 
pattern of surface temperature response from global aerosol trends 
closely matches that from warming by the well-mixed greenhouse 
gases (Harvey, 2004; Levy et al., 2008; Shindell et al., 2008; Koch 
et al., 2009; Kloster et al., 2010). A recent U.S. government report 
concluded that regional control strategies for short-lived species such 
as aerosols would have mostly global (rather than regional) impacts 
on climate (Climate Change Science Program, 2008). However, this 
conclusion does not seem consistent with the observed correlations 
between temperature and aerosol trends in source regions. Some 
GCMs do show strong regional cooling effects from anthropogenic 
aerosols, such as a 0.4 K annual mean cooling over India (Ramanathan 
et al., 2005) and a 0.5-1.0 K mean summertime cooling over China 
(Menon et al., 2002). Using a regional model, Jacobson et al. (2007) 
calculated a 0.5 K cooling over Los Angeles in winter due to anthro- 
pogenic aerosols. 

Here we use the Goddard Institute for Space Studies (GISS) GCM 3 
(Rind et al., 2007) to compare 2010-2050 temperature trends in 
a simulation with present-day U.S. aerosols vs. a simulation with no 
U.S. aerosols. We focus on only the direct aerosol effect. Our simple 
approach of zeroing AOD permits detailed analysis of a suite of 
meteorological changes, including feedback processes and the sea- 
sonality of the climate response. Also, by performing a transient 
model simulation with changing greenhouse gases, we test whether 
the interaction of aerosol and greenhouse gas effects has importance. 


2. Approach 

We use the GISS GCM 3 with 4° x 5° horizontal resolution and 
23 levels in the vertical (Rind et al., 2007). The climate sensitivity of 
this GCM to a doubling of C0 2 is 2.8 K (Rind et al., 2007), within the 
2.1 -4.4 K range of current models (Randall et al., 2007). Heat 
exchange with the ocean uses a Q-flux scheme, in which sea surface 
temperatures respond to changing climate but the horizontal 
transport of heat in the ocean is set to present-day values (Hansen 
et al., 1988). The scheme allows for diffusive mixing of heat into the 
deep ocean (Hansen et al., 2002). The Q-flux approach is widely 
used to perform ensemble and sensitivity studies (Jones et al., 
2007; Ming and Ramaswamy, 2009). The forcing from U.S. aero- 
sols is sufficiently small that it should not perturb ocean circulation 
significantly. 

For present-day aerosol loading, we use monthly mean 1990s 
AODs at 550 nm from Hansen et al. (2002); the annual mean AOD 
is shown in Fig. la. Aerosol mass concentrations over the United 
States have declined since the 1990s (Smith et al., 2011; 
Environmental Protection Agency, 2010; Leibensperger et al., 
submitted for publication), and so the climate response we calcu- 
late here may have already begun to take place. Over the United 
States, the 1990s AOD in the model averages 0.2, of which 90% is 
anthropogenic. About half the AOD is due to sulfate aerosols; the 
other main components are organic carbon (~40%) and black 
carbon ( ~ 10%). The mean single-scattering albedo is 0.93. Only the 
direct radiative effect of aerosols is included in the model. Our 
understanding of indirect effects involving aerosol— cloud interac- 
tions is poor (Forster et al., 2007) and not including them makes our 
study conservative. Philipona et al. (2009), in their analysis of 
surface radiative fluxes over Europe, found no significant trends in 
cloud cover linked to changing AOD. 

Fig. lb shows the annually averaged, top-of-atmosphere (TOA) 
direct radiative forcing from zeroing present-day aerosols over the 
United States, calculated as in Mickley et al. (1999). The forcing is 
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Fig. 1 . (a) Annual mean aerosol optical depths (AODs) at 550 nm from Hansen et al. (2002), taken as present-day conditions in the GISS GCM3 simulation, (b) Annual mean top-of- 
atmosphere (TOA) direct radiative forcing from removal of present-day aerosols over the United States, (c) 2010-2050 surface warming from increasing greenhouse gases in the 
GISS GCM 3 following the Nakicenovic et al. (2000) A1B scenario and with constant aerosols. Results are plotted as the difference between annual mean 2041-2050 and 2000-2009 
temperatures, (d) Additional 2010-2050 annual mean warming from removal of US aerosols. White areas show statistically insignificant temperature change. 
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positive, due to the removal of cooling aerosols, and is mainly in the 
shortwave. It averages 2.2 W m~ 2 over the United States, within the 
range of published studies (e.g., 2 ± 0.5 W m -2 in the multi-model 
study of Schulz et al. (2006)) and comparable to that derived from 
satellite and ground-based measurements (Chung et al., 2005). 
Over the Northeast and California the forcing is as large as 5 W nrr 2 , 
again in good agreement with observed values (Chung et al., 2005). 
The surface forcing averages 5.2 W m~ 2 over the United States, 
larger than the TOA forcing because of aerosol absorption of solar 
radiation, and is as large as 8 W m -2 over the Northeast and Cal- 
ifornia (not shown). 

We perform 2010-2050 transient-climate simulations for the 
four cases in Table 1, including calculations for 3 ensemble 
members for each case. Control Simulation 1 follows the IPCC Al B 
scenario for the well-mixed greenhouse gases, with AOD fixed 
at present-day levels. The A1B scenario is characterized by rapid 
economic growth, new energy-efficient technologies, and energy 
generation from a balanced mix of fossil fuel and alternative 
sources (Nakicenovic et al., 2000). Control Simulation 2 uses fixed 
2010 greenhouse gases and present-day AOD. Sensitivity Simu- 
lation 1 follows the A1B scenario for well-mixed greenhouse 
gases but with the AOD over the United States set to zero starting 
in 2010. Sensitivity Simulation 2 uses fixed 2010 greenhouse 
gases as in Control Simulation 2, but zeroed AOD as in Sensitivity 
Simulation 1. 

Our approach does not consider the radiative effect of U.S. 
aerosols transported outside the United States, and does not 
separate natural and anthropogenic contributions to the AOD. Most 
of the AOD from U.S. aerosol sources is over the United States itself 
(Fig. la), reflecting the short lifetime of aerosols against rainout. As 
stated above, anthropogenic sources account for 90% of the AOD 
over the United States in the Hansen et al. (2002) data. As a test, we 
conducted a 2010-2050 transient-climate simulation in which we 
zeroed only the anthropogenic component of the AOD. We found 
no significant differences in regional climate response relative to 
the case in which we zeroed total AOD. 

All simulations are initialized in 2010 with results from 
a previous transient 1950-2010 climate simulation using the same 
GCM (Wu et al., 2008), and ensemble members are generated by 
varying the initial set of random numbers used in the calculation of 
fractional cloud cover. Unless otherwise specified, results for each 
case are reported as averages of the ensemble members. 

We used the two-tailed t - test to determine the significance of 
differences between the simulations, with p < 0.05 as a threshold 
for significance. To calculate variances over the 41 model years, we 
used the residuals obtained as the difference between the time 


Table 1 

2010-2050 transient-climate GCM simulations. 



Scenario for 
well-mixed 
greenhouse 
gases 

Aerosol 
optical depths 

ATsurface, K* 

Global Eastern U.S. b 

Control Simulation 1 

A1B 

Present-day c 

1.4 

1.5 ± 0.2 

Control Simulation 2 

2010 (fixed) 

Present-day 

0.5 

0.5 ± 0.2 

Sensitivity Simulation 1 

A1B 

Zero over US, 

present-day 

elsewhere 

1.4 

2.0 ± 0.2 

Sensitivity Simulation 2 

2010 (fixed) 

Zero over US, 

present-day 

elsewhere 

0.5 

1.0 ± 0.1 


a Annual mean surface temperature change from 2000-2009 to 2041-2050. 
Standard deviations are from the three ensemble members, and are less than 0.1 K 
for global differences. 
b Region defined by 90-70W, 28-48N. 
c As shown in Fig. la. 


series of each ensemble member and the 9-year moving average for 
the ensemble. (For the first and last 4 years of the time series, 
we used the average of the first and last 9 years, respectively.) 
In calculating significance, we also adjusted the sample size, or 
degrees of freedom, to account for autocorrelation within each 
ensemble member. For this, we followed the formulation of Wilks 
(2006): 

n 

^adjusted = ^""1 1 — 0L {)/ (1 + $/), 

i' = l 

where N a dj US ted is the adjusted degrees of freedom, n = 3 is the 
number of ensemble members, N; is the sample size for each 
ensemble member (41 years), and ct; is the lag-1 year autocorrela- 
tion coefficient of the residuals of the 9-year moving average of the 
time series for each member. For surface temperature, the time 
scale of autocorrelation of the 9-year moving average is about one 
year. 

3. Climate response to removal of U.S. aerosols 

All four simulations in 2010 are out of radiative equilibrium, 
with some residual warming from the 1950-2010 spin-up stored in 
the ocean. By 2041—2050, Control Simulation 2 with constant 
greenhouse gases shows a global surface warming of 0.5 K relative 
to the spin-up years, 2000-2009. By contrast, Control Simulation 1 
with increasing greenhouse gases shows a global surface warming 
of 1.4 K by 2041-2050, as shown in Fig. lc. As is typical of future 
climate calculations, warming is strongest over high northern 
latitudes (+2-3 I<). Warming over the contiguous United States 
averages about 1 K. 

Fig. Id shows the additional surface warming caused by removing 
aerosols over the United States, as obtained by difference between the 
mean 2010—2050 surface temperatures in the Sensitivity Simulation 
1 and Control Simulation 1. This additional warming is mainly local- 
ized over the eastern United States, strongly correlated with the 
pattern of forcing (Fig. lb). It reaches 0.4-0.6 K in that region, with 
a maximum of 0.6-0.8 K in the late summer/early fall timeframe, as 
will be discussed below. About 0.2 K of significant warming occurs 
downwind over Europe and also upwind in the region of the Pacific 
High. Surface cooling or warming in isolated gridboxes elsewhere 
across the Northern Hemisphere is likely due to noise. White areas 
indicate regions where temperature changes were not statistically 
significant. In general, the pattern of surface warming is in accordance 
with the regional character of the radiative forcing. 

In Fig. 2, we present the time series of annual mean surface 
temperatures for Control Simulation 1 and Sensitivity Simulation 1 
averaged over the eastern United States. Temperatures are shown 
as the anomalies relative to the 2000-2009 mean. There is large 
interannual variability: the standard deviation of the residuals of 
the 9-year moving average for any one ensemble member is ~ 1 K. 
The bold curves in Fig. 2 show the running 9-year means of the 
ensemble averages for the control and sensitivity simulations. In 
both cases the annual mean temperatures climb by about 0.8 K 
during the 40-year simulation, due to rising greenhouse gases, but 
the simulation with zero U.S. aerosols shows a positive offset of 
0.5 K that develops immediately and persists over the course of the 
simulation. The offset is highly significant (p < 0.01 ), as determined 
by Student’s t- test. 

Fig. 3 shows the 2010-2050 annual mean response to the 
removal of aerosol for three measures of temperature over the U.S. 
domain. Fig. 3a, which shows the change of surface mean 
temperature, is the same as Fig. Id, but limited to the United States. 
We repeat the plot here to facilitate comparison with the change in 
mean daily maximum temperatures (Fig. 3b). We find that the daily 
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Fig. 2. Time series of annual mean surface temperatures anomalies over the eastern 
United States (90-70W, 28-48N) in the 2010-2050 A1B atmosphere relative to the 
2000-2009 mean. The green curves represent Control Simulation 1 with fixed 
present-day AODs; the red curves represent the Sensitivity Simulation with U.S. 
aerosols removed from 2010 on. Dotted curves show results for the 3 ensemble 
members of each simulation. Solid curves show the running 9-year means averaged 
over the ensemble members. 


maximum temperatures show a greater sensitivity to zeroed AOD 
than mean daily temperatures, with increases of 0.2-1.0 K over 
large regions in the eastern half of the United States and in the 
West, due to enhanced solar heating under clear skies. The 
response in the mean diurnal temperature range (Fig. 3c) is 
generally consistent with the change in daily maximum tempera- 
tures, with increases of 0.1 -0.2 K across much of the northern 
United States, the West, and Texas. Comparison of Fig. 3b and c 
indicates that nighttime minimum temperatures change consid- 
erably less than the mean daily maximum temperatures. 

The removal of U.S. aerosol also perturbs the hydrological cycle 
through strong surface heating and lofting of moist air. Although 
annual mean total cloud cover remains unchanged, Fig. 4a shows 
that moist convective cloud cover increases in Sensitivity Simula- 
tion 1 by 0.5-1.0% relative to the Control over much of central and 
eastern United States and nearby coastal waters. Annual mean 
precipitation also increases by 0.2-0.3 mm d 1 in a similar spatial 
pattern (Fig. 4b). The greatest enhancement in precipitation occurs 
in spring in the Northeast (0.4 mm d \ not shown). In the South- 
east, springtime precipitation increases by ~0.7 mm d -1 , but only 
during extreme events - i.e., those days with total rainfall greater 
than two sigma from the seasonal mean. 

In Fig. 5 we present the seasonality of the response to zeroed 
AOD for a suite of meteorological variables averaged over the 
northeast quadrant of the United States (box in Fig. 3a), where 
aerosol forcing is largest. This region, which we will henceforth call 
“Northeast,” covers the mid-Atlantic states south to Virginia and 
extends westward as far as Iowa and Missouri. Fig. 5a shows that 
the increase in temperature is maximum in September (0.8 K for 
daily mean and 1.0 1< for daily maximum temperatures). The change 
in the diurnal temperature range to U.S. aerosol removal is insig- 
nificant in spring but rises to 0.3 K in August-October, indicating 
large daytime heating under clear skies. In the rest of the panels of 
Fig. 5, we examine the positive feedbacks that result in the strong 
enhancement of daytime surface temperatures in late summer/ 
early fall. 

Fig. 5c shows that the increase in solar radiation absorbed at the 
surface due to zeroing of U.S. aerosols is maximum in summer at 
nearly 15 W itT 2 , as expected, but these large values persist 
through September, even as top-of-atmosphere insolation declines. 
The added solar heating enhances surface latent and sensible heat 


3 Change in mean daily temperature 



b Change in mean daily max temperature 



q Change in mean diurnal temperature range 



Fig. 3. Effects of removal of U.S. aerosols on (a) mean daily surface temperatures, (b) 
daily maximum surface temperatures, and (c) diurnal temperature range. Plots show 
the mean differences averaged over the 2010-2050 period between Sensitivity 
Simulation 1 with zeroed AOD and Control Simulation 1 with fixed present-day 
AOD. Panel (a) is the same as Fig. Id, limited to the U.S. domain. White areas 
show statistically insignificant differences. Color bars for Panels (a) and (b) are the 
same. 


fluxes by as much as 7—8 W m -2 and outgoing longwave radiation 
by 4 W nrr 2 (Fig. 5d). The change in latent heat flux peaks sharply in 
June, while the responses in sensible heat and longwave fluxes are 
greatest three months later in September, when declining levels of 
soil moisture diminish the importance of the latent heat flux in the 
surface energy budget. While such a loss of soil moisture occurs 
normally in late summer/early fall even in the absence of forcing, 
the enhanced solar heating in the simulation with zeroed AOD 
further reduces the available soil moisture in this timeframe, by 
about 4% (Fig. 5e). The depletion in soil moisture occurs despite the 
annual mean increase in precipitation rates (Fig. 4b), most of which 
occurs in spring. This depletion likely shifts an even larger portion 
of the outgoing surface energy toward sensible heat in late summer 
and early fall. The enhancement in sensible heat more strongly 
affects local surface temperatures than an equivalent increment in 
latent heat flux, since the latter can efficiently transport the added 
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Fig. 4. Effect of removing U.S. aerosols on annual mean (a) moist convective cloud cover and (b) precipitation rates. The plots represent the mean 2010-2050 differences between 
Sensitivity Simulation 1 with US aerosols removed and Control Simulation 1 with fixed present-day AOD. White areas indicate statistically insignificant differences in precipitation. 


heat elsewhere (Boer, 1993; Brabson et al., 2005). Thus the decline 
in soil moisture contributes to a positive feedback on surface 
temperatures in late summer/early fall in the Northeast. By 
contrast, the Southeast United States does not experience signifi- 
cant declines in soil moisture availability. 

This positive feedback in the Northeast is further augmented by 
changes in surface relative humidity and low cloud cover (Fig. 5f). 
Relative humidity at the surface exhibits a large 2% decline in the 
July-October timeframe owing to warmer temperatures and drier 
soils when U.S. aerosols are removed (Dai et al., 1999). The trend in 
low cloud cover follows that of surface relative humidity, as 
expected from observations during the warm season over the 
continental United States (Del Genio and Wolf, 2000). The large 
decline in low cloud cover in late summer/early fall in Sensitivity 
simulation 1 allows an even greater flux of solar energy to reach the 
surface, as can be seen in Fig. 5c. The net effect is to sustain the 
increase in surface temperatures and the diurnal temperature 
range through early fall (Fig. 5a and b), even as the direct forcing of 
aerosol declines. 

In Fig. 6, we further examine the regional meteorological 
response to zeroed AOD in the July— September timeframe, when 
surface temperatures in the Northeast are most sensitive to the 
change in aerosol loading. Fig. 6a shows that mean July-September 
sea level pressures across the East drop by 2-4 hPa, owing to added 
surface heating. The change in sea level pressure triggers increased 
cyclonic circulation in the 900 hPa wind fields. The pattern of 
response in moisture transport follows that of the 900 hPa winds, 
with increased offshore transport of moisture from the Northeast. 
Although precipitation in Sensitivity Simulation 1 shows strong 
enhancement across the Northeast when averaged annually 
(Fig. 4b), in the July-September timeframe mean rainfall rates 
increase significantly only over nearby coastal waters, further 
drying out the Northeast (Fig. 6b). 

We next investigate the fate of the additional heat generated at 
the Earth’s surface, again in the July-September timeframe. Fig. 7 
shows the profiles of the mean July-September 2041-2050 
temperature changes over the Northeast United States for the 
Control Simulation 1 with fixed aerosols and the Sensitivity Simu- 
lation 1 with zeroed U.S. aerosol. The changes are calculated with 
respect to the 2000-2009 mean of Control Simulation 1. Warming 
from the A1B trend in greenhouse gases peaks at 2.0 K in the mid- 
troposphere ( ~380 hPa), as is typical in other model studies (Meehl 
et al., 2007). The profile of additional heating caused by zeroed AOD 
shows a steep vertical gradient in the lower troposphere, with 
a temperature enhancement of 0.6 K at the surface, compared to 
- 0.2 K at 660 hPa. The overlying stratosphere cools by ~ 0.4 K in the 
simulation with zeroed AOD relative to the Control during this time 
period. The cooling is likely caused by a ~2% enhancement in 


specific humidity in the mid-troposphere (p < 0.10, not shown) and 
the concomitant increase in infrared opacity. 

Fig. 8 shows the mean July-September temperature response to 
zeroed U.S. aerosols over a larger domain at 500 hPa. At these 
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Fig. 5. Seasonal dependence in the response of several variables to removal of U.S. 
aerosol in the A1B atmosphere, averaged over the Northeast (box in Fig. 3a). Shown are 
the monthly mean differences averaged over the 2010-2050 period between Sensi- 
tivity Simulation 1 and Control Simulation 1 for (a) daily mean and daily maximum 
surface temperatures; (b) diurnal range in surface temperature; (c) absorbed solar flux 
at the surface; (d) surface fluxes of latent heat, sensible heat, and longwave radiation; 
(e) change in soil moisture availability; (f) low cloud cover and surface relative 
humidity. Positive values in Panel (d) denote increased upward flux. Percentage values 
in Panel (f) denote changes in absolute amount. Error bars show the standard deviation 
across the 41 model years and the 3 ensemble members for each case. Error bars for 
different curves are offset for clarity. 
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Fig. 6. Effect of removing U.S. aerosols on mean July-August-September (a) sea level pressure (color background) and 900 hPa winds (arrows) and (b) precipitation rates 
(color background) and vertically integrated moisture transport (arrows). The plots represent the differences between Sensitivity Simulation 1 with U.S. aerosols removed and 
Control Simulation 1 with fixed present-day AOD. The change in moisture transport is summed through the depth of the troposphere. White areas indicate statistically 
insignificant differences in sea level pressure and precipitation. To show the patterns in the changes in 900 hPa winds and moisture transport, we plot all values regardless of 
significance. 


altitudes, the prevailing westerlies carry the warm air in a band 
stretching across the north Atlantic to Europe, with temperatures in 
the band averaging about 0.2-0.3 K greater than in the Control. 
Over Europe the warm air subsides, heating the surface there by 
about 0.5 K in the July-September timeframe, or about 0.2 I< in the 
annual mean (Fig. Id). Fig. 8 also shows a northward shift in the 
westerlies, as heating over the United States steepens the meridi- 
onal temperature gradient poleward, but flattens this gradient 
toward the south. 

We find that the removal of U.S. aerosol intensifies summer heat 
waves in the Northeast. Fig. 9 shows the 2010-2050 frequency 
distributions of daily mean surface temperature anomalies over the 
Northeast for July and August, when heat waves are most common. 
(Daily maximum temperatures are available only as monthly 
averages.) In Sensitivity Simulation 1 (red), with zero AOD over the 
United States, the median temperature increases by about 0.5 K 
relative to Control Simulation 1 (green), consistent with the change 
in mean July-August temperatures. But temperatures during heat 
waves, at the high end of the distribution, are 1—2 K higher in the 
simulation with aerosol removed. This effect is not seen in winter 
when aerosol radiative forcing is minimal. 

In order to test whether the aerosol warming effect is contingent 
on Al B growth in greenhouse gases, we compared Control Simulation 



Temperature (K) 

Fig. 7. Vertical profiles of the mean response of temperature to changes in greenhouse 
gases and aerosols, averaged over the Northeast United States (boxed region in Fig. 3a) 
for Jul-Aug-Sept 2041-2050, relative to 2000-2009. The bold curve represents the 
effect of increasing Al B greenhouse gases with constant aerosols (Control Simulation 1 ). 
The dashed curve represents the effect of the same Al B increase in greenhouse gases but 
with zeroed US aerosol (Sensitivity Simulation 1 ). Error bars show the standard devia- 
tion across the 41 model years and the 3 ensemble members for each case. Error bars 
have been offset for clarity. Circles denote an insignificant difference between the 
temperature responses in the two simulations. 


2 (constant greenhouse gases and present-day AOD) with Sensitivity 
Simulation 2 (constant greenhouse gases but zeroed AOD over the 
United States). Both these simulations reach radiative equilibrium by 
- 2025. We find that the climate response to aerosol removal in the 
test simulation is similar to the response we calculate with changing 
greenhouse gases (Table 1). As before, the surface temperature 
response is immediate and an offset of 0.5 K remains throughout the 
2010-2050 time period. 

4. Discussion 

Previous research examining the regional climate response to 
aerosol forcing has mainly focused on Asia, where the large 
aerosol burden has been estimated to reduce solar insolation by 
-15-25 W rrT 2 in some regions (Krishnan and Ramanathan, 2002; 
Ramana et al., 2004). Although we calculate a much weaker surface 
aerosol forcing over the United States in our model, about -5 W m -2 
as an annual mean, we find that removal of U.S. aerosol leads to 
significant regional warming, with annual mean surface temperatures 
increasing by 0.4— 0.6 K in the East. During summertime heat waves in 
the Northeast, mean daily temperatures rise by as much as 1-2 K. 
In late summer-early fall, the diurnal temperature range increases 



Fig. 8. Effect of removing U.S. aerosols on mean July-Aug-Sept 500-hPa temper- 
atures (color background) and wind fields (arrows). The plots represent the 
mean 2010-2050 differences between Sensitivity Simulation 1 with U.S. aerosols 
removed and Control Simulation 1 with fixed present-day AOD. White areas indi- 
cate statistically insignificant differences in temperature. For clarity, all changes in 
windspeed greater than 0.3 m s -1 are shown, regardless of the significance of the 
difference. 
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Fig. 9. Cumulative probability distributions of daily mean surface temperatures 
anomalies over the Northeast during July- August for the 2041-2050 period. Results 
are shown for Control Simulation 1 (green), with changing A1B greenhouse gases and 
fixed present-day aerosol burden, and for Sensitivity Simulation 1 (red), with changing 
A1B greenhouse gases but zeroed aerosol over the United States. Each point represents 
the daily mean temperature anomaly relative to the median of the Control, averaged 
over the Northeast for the 9 simulation years and the 3 ensemble members. The x-axis 
is a probability scale such that a normal probability distribution would plot as 
a straight line. 

by 0.3 K in this region, consistent with strong daytime heating 
amplified by the loss of low cloud cover and depleted levels of soil 
moisture. 

Our model results suggest that recent declines in U.S. aerosol 
sources (Streets et al., 2004; Smith et al., 2011) may have contrib- 
uted to the observed 0.2-0.5 K warming over the United States in 
the last decade (2001-2010), relative to the previous 20 years 
[http://data.giss.nasa.gov]. Our finding of an increase in the diurnal 
temperature range when aerosol is removed is consistent with 
Makowski et al. (2008), who reported anticorrelations in annual 
mean SO 2 emissions and the observed diurnal temperature ranges 
for several sites in Europe. For the United States, observations show 
a general increase in the diurnal temperature range since the 1950s 
(e.g., Vose et al., 2005; Zhou et al., 2009), but to our knowledge no 
study has examined decadal relationships in this variable and U.S. 
aerosol loading. Evidence of weekly cycles in the observed diurnal 
temperature range over many sites in the United States (Forster and 
Solomon, 2002) imply that anthropogenic emissions can influence 
daily maximum or minimum temperatures. 

We identify a positive feedback occurring in the Northeast in 
late summer/early fall that may explain the large sensitivity of 
regional temperatures to aerosol forcing in our results. In this 
feedback, strong surface heating depletes soil moisture by 
midsummer, which enhances the outgoing surface fluxes of 
sensible heat and longwave radiation through early fall. A similar 
soil moisture feedback has been implicated in the 2003 European 
heatwave (Schar et al., 2004). Warmer temperatures and drier soils 
also reduce low cloud cover in the July— September timeframe, 
allowing more sunlight to reach the surface and increasing 
temperatures still further (Dai et al., 1999; Dai et al., 2006). These 
feedback mechanisms illustrate the importance of soil processes in 
governing regional climate response. 

The late summer/early fall loss of soil moisture occurs with 
zeroed AOD despite an increase in annual mean precipitation of 
~0.2-0.3 mm d -1 across much of the United States. In global 
projections of 21st century rainfall, the United States falls in 
a transition zone, with most models agreeing that precipitation will 
increase over latitudes north of the United States and will decrease 
to the south (Christensen et al., 2007). In the transition zone, 
models yield a more mixed picture. Our result suggests that 
declining aerosol sources over industrial regions in the mid- 
latitudes may tilt the balance toward increased rainfall. 


We find that the spatial pattern of surface warming over the 
United States strongly corresponds to the pattern of aerosol forcing. 
While this result is consistent with observed relationships between 
forcing and temperature response (Ruckstuhl et al., 2008; Krishnan 
and Ramanathan, 2002), it differs with some model studies report- 
ing little spatial correlation between regional radiative forcing and 
surface temperature response (Brasseur and Roeckner, 2005; 
Shindell et al., 2007; Levy et al., 2008; Shindell et al., 2008; Liao 
et al., 2009). These model studies however did not separate the 
effects of aerosols and tropospheric ozone, combining them as an 
ensemble of short-lived climate forcing agents. Temperature 
response to ozone forcing is mainly on the hemispheric scale 
(Mickley et al., 2004; Yoshimori and Broccoli, 2008). In Asia and other 
regions of the developing world where black carbon aerosol is 
relatively abundant, changes in aerosol loadings lead not just to local 
cooling (Krishnan and Ramanathan, 2002), but also to diffuse 
warming since black carbon, like ozone, heats the mid- to upper 
troposphere (Menon et al., 2002; Zhang et al., 2009; Jones et al., 
2011 ). By contrast, black carbon over the United States in our study 
accounted for only about 10% of the AOD. The indirect effect of 
aerosols on clouds, which was not included in our work, could 
further diffuse the climate response. In addition, the soil moisture 
feedback, which in our study intensified the local temperature 
response in the Northeast, may not function as effectively elsewhere. 
In regions with ample soil moisture, increased insolation would have 
less impact on temperatures (Clark et al., 2010). 

A shortcoming of our approach is that we did not take into 
account the transport of U.S. aerosol elsewhere. In their model 
study, Fischer-Bruns et al. (2010) reported that Arctic surface 
temperatures are sensitive to changes in North American aerosol 
sources, and some of the previously mentioned studies also found 
a strong Arctic response to global trends in short-lived species (e.g., 
Levy et al., 2008; Shindell et al., 2008). However, Arctic aerosol 
pollution is dominated by transport from Asia and Europe, with U.S. 
sources having less influence (Stohl, 2006; Hirdman et al., 2010). 
Polluted air masses leaving the east coast of the United States 
experience lofting and diabatic heating on the warm conveyor belts 
prevalent in the North Atlantic, a process which washes out most 
aerosols (Law and Stohl, 2007; Fisher et al., 2010). The large Arctic 
temperature response calculated by Fischer-Bruns et al. (2010) may 
be partly due to the effect of global volcanic emissions, which 
together with North American anthropogenic aerosol sources are 
set to zero in their sensitivity simulation. Using a similar model as 
Fischer-Bruns et al. (2010), Graf et al. (1997) found that volcanic 
aerosol imposes forcings as great as 3 W nrT 2 over high northern 
latitudes. 

Another shortcoming of this study is that the coarse model 
resolution we used cannot resolve orographic effects on meteo- 
rology. These effects would have greater importance in the western 
United States, however, than in the East. Our study also did not 
include a realistic variation of aerosol sources over time. Though 
simple, our approach of abruptly zeroing U.S. aerosol allows for 
greater clarity in understanding the suite of meteorological changes 
and feedbacks. A follow-up study (Leibensperger et al., submitted 
for publication) uses the same climate model to calculate the 
radiative forcing and climate response to historical (1950-2000) 
and projected (2000-2050) trends in anthropogenic aerosol 
sources. 

The response to changing aerosol sources is likely to be model- 
dependent, and the robustness of our results needs to be tested 
with other GCMs. There is a clear need for concerted multi-model 
and observational studies to investigate relationships between 
regional aerosol radiative forcing and climate response. In partic- 
ular, the sensitivity of key feedback variables such as soil moisture 
and low cloud cover to changing conditions needs to be better 
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understood (Clark et al., 2010; Williams et al., 2006; Soden and 
Held, 2006). With these caveats, our results suggest that U.S. 
efforts to improve air quality and limit acid rain through regulation 
of aerosol sources could significantly exacerbate regional climate 
warming. Our work underscores the importance of integrating air 
quality and climate concerns in the design of emission control 
strategies. 
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